Highly anisotropic Bose-Einstein condensates: crossover to lower dimensionality 
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We develop a simple analytical model based on a variational method to explain the properties of 
trapped cylindrically symmetric Bose-Einstein condensates (BEC) of varying degrees of anisotropy 
well into regimes of effective one dimension (ID) and effective two dimension (2D). Our results are 
accurate in regimes where the Thomas-Fermi approximation breaks down and they are shown to be 
in agreement with recent experimental data. 
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I. INTRODUCTION 

There has been a growing interest in Bose Einstein 
condensation in effective lower dimensions, with exper- 
imental realization finally becoming reality M. Apart 
from the purely academic interest, there is a very prac- 
tical interest as well arising from the rapidly developing 
endeavor to create efficient atom waveguides with poten- 
tial applications in interferometry and gyroscopes. For 
trapped gases effective lower dimensionality means that 
excitations along the tightly confined dimcnsion(s) are 
energetically not allowed. This is a limiting case of highly 
anisotropic condensates which are becoming more com- 
mon with a new generation of BEC experiments on sur- 
face micro traps |3, 0] . 

Interest in condensates of extreme anisotropy is evident 
in numerous current theoretical and experimental work 
on atomic waveguides, quasi-2D configurations [0, O and 
the Tonks-Girardeau (TG) ID limit g |, § in which im- 
penetrable bosons show fermionic properties. Magnetic 
waveguides for neutral atoms of diverse design have been 
constructed in many laboratories J3|, [|, H], [fl], O, [J3| . 
Bessel beams |14J have been used to produce wave-guide- 
like optical confinement. With a slight longitudinal po- 
tential these waveguides can be treated as high aspect 
ratio traps. Atoms have been trapped in a 2D optical 
lattice consisting of quasi- ID optical wells |L5| each with 
an aspect ratio of up to about 2000. Work is under- 
way at JILA on various waveguide and near-waveguide 
configurations fll6f . On the theoretical side the various 
regimes of quantum degeneracy in both the ID and 2D 
limits have been studied in detail in several recent papers 

But a simple model is lacking which would describe 
how a condensate changes as it becomes more anisotropic 
and eventually crosses over to effective lower dimension- 
ality. For condensates in 3D, of atoms with positive scat- 
tering length, a highly effective and yet essentially simple 
analytic description was achieved through the Thomas- 
Fermi Approximation (TEA) (p7|. It neglects the kinetic 



energy in the Gross-Pitaevskii (GP) equation which is 
the mean-field description of the condensate and it is 
justified for large condensates with aspect ratios of the 
order of unity. But as the aspect ratio deviates farther 
from unity, the kinetic energy in the constricted direction 
becomes increasingly more important and the TFA does 
not work so well. Although there have been numerous, 
insightful refinements of the TFA ||, (|| |o[ |l]> || they 
do not directly improve on it for highly anisotropic BEC. 
It is desirable to have a theoretical model comparable to 
the TFA in simplicity but one which is successful in de- 
scribing condensates from the 3D regime with increasing 
degree of anisotropy all the way to regimes of effective 
lower dimensionality. That is the objective of this paper. 
For longitudinally homogeneous systems, we have re- 
cently coauthored a study of the crossover from 3D to 
effective ID [f53|. Here we will consider gases harmoni- 
cally confined in all directions and a crossover to effective 
2D as well. We will use a variational approach to obtain 
analytic expressions for the chemical potential, total en- 
ergy, release energy and also the energy spectrum valid 
over a wide range of parameters that include the regime 
of crossover to lower dimensions. In Sec. II we compare 
our chemical potential with accurate numerical solutions 
of the GP equat ion and the Thomas-Fermi expression 
and then in Sec. Ill we discuss the crossover regime. In 
Sec. IV we compare the release energy from our model 
with experimental data from Ref. |jj. Then in Sec. [V] we 
consider the Bogoliubov equations for quasiparticle ex- 
citations and obtain the energy spectrum for geometries 
close to effective ID and 2D. 



II. CHEMICAL POTENTIAL 

At zero temperature in the lowest order mean-field 
approximation the condensate is described by the wave 
function that minimizes the GP energy functional |27[ ] 
which in cylindrical co-ordinates is 
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We have defined the aspect ratio as the ratio of the ra- 
dial and the axial trapping frequencies 7 = uj r /uj z . The 
interaction strength is determined by the scaled scatter- 
ing length g = a/ao and by the atom number N which 
we take to be equal to the condensate number assuming 
temperatures close to absolute zero. The co-ordinates 
r and z are scaled by the respective oscillator lengths 
ro = y^h/mujr and zq — ^J%/mu: Z) with the mean oscil- 
lator length a = y/fi/muj being defined by the geometric 
mean frequency, uj = (wjw z )'' 3 . The condensate wave- 
function is (e~"" e /v / 27r) ( f > (?', z) with m = and normal- 
ization J rdr J dz\&(r, z)\ 2 = 1. 

Functional minimization 6E/5Q = with the normal- 
ization condition gives the GP equation 
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This determines the optimum condensate wavefunction 
$ and the corresponding chemical potential fi. For con- 
densates with all its spatial dimensions of comparable 
magnitudes the Thomas-Fermi approximation provides a 
accurate value of the chemical potential 



, TF = *£ (15Ngf 5 . 
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This expression has no dependence on the aspect ra- 
tio, whereas we would expect the chemical potential to 
change as the aspect ratio changes. In order to obtain 
analytic expressions for $ and /j, that will have the cor- 
rect dependence on the aspect ratio we take a variational 
approach with a trial wavefunction containing a few pa- 
rameters and then minimize the energy functional with 
respect to them [p3. In choosing our trial function we 
note that for highly anisotropic traps, in the direction 
of weak confinement the condensate size far exceeds the 
oscillator length and the kinetic energy (which scales as 
the inverse square of the size) becomes negligible, so a 
Thomas-Fermi form is appropriate. In the squeezed coor- 
dinate the interaction term has lesser relative importance 
so that a Gaussian form is suitable. 



A. Cigar Geometry 

For the cigar geometry with 7 > 1, we take our nor- 
malized trial function to be 
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where f3 r and d are variational parameters and Q(x) is 
the unit step function. In evaluating the energy func- 
tional, Eq. (jj), for this trial function, we recognize that 
the kinetic energy in the transverse direction should be 
included, but not that in the axial direction. We thus 
obtain 
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FIG. 1: Comparison of chemical potential obtained from var- 
ious approaches, for cigar-shaped traps with (a) 7 = 100 and 
(b) 7 =1000 . The methods are numerical solution of the 
Gross-Pitaevskii equation: the Thomas-Fermi approximation; 
analytic expression Eq. (|7J) with optimized f} r and d\ and the 
same with /3 r — 1. The optimum value of /3 r is plotted with 
the right axis scale. 



Minimizing E with respect to d and r gives 
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These coupled equations can be easily solved numerically, 
and the chemical potential corresponding to the optimum 
parameters is 



fj, = huj 
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This satisfies the thermodynamic relation /i = dE/dN 
as expected. It is apparent from the unsealed form that 

2/3 

the last term is oc u> z and hence vanishes in the limit 
lu z — > 0, also in that limit r — > 1 so that the chemical 
potential fj, equals the transverse ground-state energy Tiuj r 
as we would expect. In fact for very elongated traps, to 
a very good approximation we can set [3 r — 1 and Eq. 
(R) reduces to a self-contained analytic expressions for fi 
dependent only on N, a and 7 all of which are measurable 
parameters of the system. 



In order to test the accuracy and validity of our ap- 
proximation we compare the analytic expressions for the 
chemical potential with accurate numerical solutions of 
the GP equation, obtained using a discrete variable repre- 
sentation (DVR) mesh in r (Laguerre DVR) and z (Her- 
mite DVR) |35], [3(| |37]]. We checked for convergence of 
the numerical results as a function of mesh size and range. 
Typically, 1,000 to 1,500 mesh points sufficed. 

In Fig. |l| we plot the chemical potential computed in 
various ways for cigar traps of two different aspect ratios. 
It is clear that the variational chemical potential using ex- 
pression (nn with optimized /3 r closely follows the numer- 
ically computed chemical potential over a large range of 
Ng. For low linear densities, the expression with /3 r = 1 
is sufficient but it fails as the density increases. On the 
other hand the Thomas-Fermi expression (0) is accurate 
at high densities and breaks down at low densities since 
it approaches zero while the correct [i should approach 
the zero-point energy. 



B. Pancake Geometry 

On carrying through a similar analysis for the pancake 
geometry with 7 <C 1, taking as trial function 
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we find that the energy functional, neglecting the trans- 
verse kinetic energy, is 
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The equations for the optimum parameters that minimize 
the energy are 
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and the corresponding expression for the chemical poten- 
tial 



fl — TlLO 



7" 



-2/3 



Pz 



1 



8/3, 



1/4 



yjNg-y 1 / 3 



(11) 



As we did for the cigar geometry, we plot the chemical po- 
tential evaluated in various ways in Fig. ||. The deviation 
of the Thomas-Fermi expression from the correct expres- 
sion for \x is more pronounced in this geometry while our 
variational expression from Eq. (|ll| ) reproduces almost 
exactly the numerical solution of the GP equation over 
the range of Ng shown. Even the analytic expression for 
fi with f3 z = 1 is accurate to large values of Ng. The rea- 
son for the better agreement is larger relative importance 
of the kinetic energy in the tightly confined direction. 
Note also that in both Fig. |] and Fig. 0, differences be- 
tween the TEA and the numerical results persist to larger 
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FIG. 2: The chemical potential for pancake shaped traps ob- 
tained from various approaches, for aspect ratios (a) 7 = 0.01 
(top) and (b)7 = 0.001. In both cases, the result with opti- 
mized b and /3 Z is essentially indistinguishable from the nu- 
merical solution of the GP equation. 



values of Ng as the anisotropy increases i.e. 7 deviates 
more from unity. 

In Fig. we plot the chemical potential obtained from 
our variational expressions in Eqs. (0) and (11) as a func- 
tion of the aspect ratio 7, over several orders of magni- 
tudes that vary from very elongated cigar condensates 
to extremely oblate pancake shapes. Strong variation 
with the aspect ratio is apparent. In the same figure 
the Thomas-Fermi chemical potential htf (0) appears 
as flat horizontal lines tangential at 7 = 1 to the lines 
corresponding to our calculated /1. Thus our expressions 
for the chemical potential agree with the Thomas-Fermi 
value for spherical symmetry. 



III. CROSSOVER TO LOWER DIMENSIONS 

The variational parameters we have used to obtain the 
chemical potential also provide a measure of the effec- 
tive dimensionality of the condensate. A trapped Bosc 
gas is considered to be in effective lower dimension if 
excitations in the tightly confined dimension are frozen. 
Thus a criteria for crossover to lower dimensionality M 
is when the interaction energy per particle is comparable 
to the energy (~ fuvught) to excite the first excited mode 
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FIG. 3: The variationally optimized chemical potential \i is 
plotted as a function of the aspect ratio 7 for different mean 
field strengths. For 7 > 1 Eq. ([]) is used and for 7 < 1 
Eq. (|ll|). The horizontal dashed lines, that represent the 
Thomas- Fermi chemical potential /j.tf, are tangential to the 
corresponding variational \i at their lowest points at 7 = 1. 



N < 0.38 x (ro/a)j. The impenetrable regime requires a 
much more stringent condition kpr^/a <C 1 |l7| , where 
kp is the Fermi wavevector of the equivalent Fermi sys- 
tem. For a large enough number of atoms and at zero 
temperature, the asymptotic form of the harmonic oscil- 
lator modes provides a estimate for the Fermi wavevector 
kp = v2N/zq, so that the condition for an impenetra- 
ble gas translates to TV <C ^j(a/ro) 2 - For typical traps 
a<ro which means that the particle number needs to be 
much smaller than that required for one-dimensionality. 
Thus it follows from our considerations above that the 
transverse profile in the impenetrable regime has to be 
that of the ground state of the transverse trap potential. 
Our variational functions however cannot be applied 
to the Tonks- Girardeau regime, since the Thomas-Fermi 
profile in the axial direction has to be replaced by a 
square-root of a parabola |18| , and the axial energy is 
simply the Fermi energy for N particles in ID. Thus while 
the energies and chemical potential we have derived can 
be applied over a wide range of parameters spanning ef- 
fective ID and 2D, they are not meant to describe the 
impenetrable regime. 



in the tightly confined direction. For a weakly interact- 
ing Bose gas at low densities the interaction energy is 
roughly equal to the Thomas-Fermi or 3D chemical po- 
tential so that the condition for lower dimensionality is 
fiTF < fraught which in our units become 
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In Figs. |] and g, we have plotted the optimum value of 
(3 r and f3 z respectively along the right axis. The values 
corresponding to the crossover to one dimension gNm = 
O.387 5 / 6 and to two dimension gN 2 D < O.387" 5 / 3 are 
indicated. We see that this roughly corresponds to (3 r ~ 
0.8 and j3 z ~ 0.8. Thus the values of j3 r and j3 z give 
a measure of the dimensionality of the system: when 
/3 r ~ 1 the system is effectively one-dimensional since 
the transverse profile of the condensate coincides with 
that of the transverse ground state, likewise when j3 z ~ 
1 the system is effectively two dimensional. As these 
parameters deviate from unity the system approaches 3D. 
It is worth noting at this point that in Ref. p3f a similar 
criterion was used to study the 1D-3D crossover for a 
homogeneous system, although the parameter used was 
different, being the effective ID interaction strength. 

The ID case is of particular interest because at extreme 
low density and tight confinement the gas becomes an im- 
penetrable Tonks-Girardeau gas [|| m]. That regime is a 
subset of the regime of one-dimensionality. A wavefunc- 
tion for a gas of such impenetrable bosons can be mapped 
to that of a gas of free fermions. 



With slight rearrangement the criterion (12) we have 
used above for one dimensionality can be written as 



IV. COMPARISON WITH EXPERIMENT 

In the experiment in Ref. M where BEC in effective 
lower dimensions was achieved, the crossover from a 3D 
to lower dimensionality was deduced (i) by observing a 
sudden change in the aspect ratio of the released con- 
densate when the number of atoms was lowered below a 
certain value, and (ii) by observing a saturation of the re- 
lease energy at the zero-point kinetic energy in the tightly 
confined direction. Since the measurement of the aspect 
ratio was done several milliseconds after release from the 
trap, a proper theoretical description would involve the 
dynamics pq ] of the expansion of the condensate which 
we do not consider in this paper. However due to conser- 
vation of energy we can estimate the release energy in our 
model and compare it with the experimentally observed 
values. 

The release energy is the energy of the system after the 
traps are switched off, which is just the sum of the kinetic 
and the interaction energies of the condensate before re- 
lease. For the cigar geometry 7 3> 1, the variational 
calculation in Sec. II gives the release energy per particle 
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where we use the optimized parameters (3 r and d from 
Eq. (||). The release energy in M is plotted as a func- 
tion of the half-length (Z) of the condensate. It turns out 
that the axial expansion in the ID experiment was negli- 
gible within the time of flight till measurement so that for 
comparison we simply need the initial half-length before 
release. 

Considering the form of our trial function fl) it is 
tempting to identify d as the half-length of the con- 



densate. However, this would not be correct since the 
trial function assumes a uniform cylindrical cross-section 
along the length of the condensate, whereas as we ap- 
proach the 3D regime, the condensate tends to be an 
ellipsoid, so that the center will tend to bulge out more 
than the ends. 

One could invoke arguments based on geometry to ex- 
tract the half length from our expression for d. However 
a much easier way to do that is to make sure that our 
expression for the half-length has the correct form in the 
ID limit and the 3D (Thomas- Fermi) limit 
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The ID limit is obtained trivially d((3 r — 1) — > Z\u- In 
considering the 3D limit we note that 1/^/37 is the width 
of transverse variational profile, and the second term in 
the expression (0) for (3 r measures the deviation from ID; 
as the condensate moves away from the ID regime and 
approaches a ellipsoidal shape, that deviation would be 
maximum at the center and negligible at the end. Thus 
we take an average value for the deviation and define the 
modified parameter f3 r 
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It is easy to see that in the 3D limit J3 r <C 1, we have 
/3 2 ~ Sdy 1 ' 3 /3Ng and we get the correct limit d(/3 r ) — > 
Z 3 £>, while we get the correct ID limit as well d(J3 r = 
1) — > Z\y). Thus in the intermediate regime we expect 
the half-length to be given to a good approximation by 



Z = 3N 9 p rl 2 / 3 . 
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We stress that the release energy should be evaluated 
with the optimized variational parameters from Eq. (|6J). 
The variational process optimizes energy and not con- 
densate dimensions and hence there is no inconsistency. 
In Fig. H we plot the release energy per particle from 
Eq. (|l3|) versus the half length from Eq. (|l6|) for the pa- 
rameters corresponding to the ID experiment in Ref. p| : 
sodium atoms in a magnetic trap with radial frequency 
uj r = 2ir x 360 Hz and axial frequency lo z = 2n x 3.5 Hz, 
so that the aspect ratio is 7 ~ 103. A comparison with 
plot (3b) presented in Ref. jjj shows that our expression 
for the release energy closely follows their measured data 
for the release energy; as the system approaches effec- 
tive ID the saturation of the release energy at the radial 
zero-point energy is clear. 

In the pancake geometry the release energy per particle 
from our variational calculations is 
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Considering the agreement in the cigar case we expect 
this expression to agree well with experimentally mea- 
sured release energy in this geometry. We will not do an 
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FIG. 4: The release energy per particle E re i/h plotted as a 
function of the half-length of a cigar-shaped condensate of 
trapping frequencies u) r — 2tt x 360 Hz and lo z = 2n x 3.5 
Hz. Compare with Fig. 3(b) in Ref. |l|. The horizonal line 
represents the transverse zero-point energy. 



actual comparison here for two reasons: First the above 
expression is for strictly cylindrical geometry whereas the 
trap configuration in the 2D experiment in [|l| was not. 
Secondly unlike the ID experiment there is significant ex- 
pansion of the condensate in all directions till the time of 
measurement which involves the dynamics of expansion 
which is not a topic of this paper. 



V. EXCITED STATES 

In Sec. II we found the chemical potential to be in- 
creasingly different from the Thomas-Fermi value as the 
anisotropy of the system increases; since the quasiparti- 
cle spectrum depends on the chemical potential we should 
expect some consequent changes. The energy spectra for 
purely ID and 2D system have been known for some- 
time p2fl , also solutions have been found for cylindrical 
geometry in 3D using TEA in all directions |59]; [hj, p~L| . 
We will consider here how the spectrum changes as the 
system deviates from effective ID or 2D. 

Consider the Bogoliubov equations for the quasiparti- 
cle excitations; as with the ground state we assume fac- 
torized spatial dependence of the quasi-particle modes 
U(r,z) and V(r,z). We first discuss the pancake geom- 
etry. The strong axial confinement makes it reasonable 
to assume that U, V have the same z-dependence as the 
condensate: 



U nm (r,z) = R+ m (r)((z), V nm (r,z) = R nm (r)((z), 
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This allows us to integrate out the axial dependence. On 
using our expression Eq. (nil) for the chemical potential 



the terms reflecting axial kinetic and potential energy 
drop out both in the GP and the Bogoliubov equations 
and they take the form 
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Here g 2 n = 6 4 /16 defines an effective transverse two- 
body interaction strength and fuv r e nm are the transverse 
normal mode energies. 

The sum and difference of the Bogoliubov equations 
can be converted to normal mode equations for the trans- 
verse density and phase fluctuations 
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This is achieved by writing the small amplitude fluctua- 
tions of the condensate wavefunction as 5(f) — ^/pe lSe — 
y/po about its equilibrium <f>$ = */po, whereby we get the 
well known forms for the linearized density and phase 
fluctuations E2Tl 
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Separating the equations and neglecting higher than sec- 
ond order derivatives, being small due to the large trans- 
verse size, we find that the phase fluctuations obey the 
same hydrodynamic equation as density fluctuations 
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where p nm {r) = 6 nm {r) = r^W nm (r). Taking the 
Thomas-Fermi expression for the transverse condensate 
density this can be reduced to a hypergeometric equation 
by defining x — r 2 /b 2 : 
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For a similar analysis of the cigar geometry, we define 

U n (r,z) = t/>{r)Z+(z), V n (r,z) = ^(r)Z"(*), 
$(r,z) = ip(r)<p(z) with ip(r) = ^Wr?-'^ l2 (26) 

and in the limit of extended axial dimension we find that 
the transverse density and phase fluctuations both satisfy 
the hydrodynamic equation 
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with p n (z) = 9 n (z) = P n (z/d) and an effective ID in- 
teraction strength g\u = d 3 /3. For the Thomas-Fermi 
condensate density in the longitudinal direction, P n (z/d) 
are Legendre polynomials with eigenvalues 
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The expressions for the normal mode energies and 
eigenstates in (|2|) and (|2S|) are identical with results 
for effective ID and 2D obtained in previous works 
However the total energy for each mode 
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has to include the energy due to the tightly confined di- 
rection, which in effective 2D is simply fiu z /2 and in ID 
is hu r /2. However as the system deviates from lower 
dimensionality we see that contribution will change and 
the total energy will be 
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with n,m — 0, 1,2, ■••. These expressions for the en- 
ergy spectrum would apply when the system has devi- 
ated form effective lower dimensionality but has not com- 
pletely attained Thomas- Fermi behavior in all directions. 
In this context we may point out that for a 3D cylindri- 
cally symmetric system well described by Thomas-Fermi 
in all dimensions, the energy spectrum is quite different 
from the ID spectrum J40|, |l|. 

Moreover the behavior of our eigenfunctions in the 
weakly confined dimension(s) depend on the behavior 
in the tightly confined dimcnsion(s) through the inter- 
dependence of the variational parameters. Thus in the 
crossover regime where the (i r and (i z start to deviate sig- 
nificantly from unity, the shapes of the eigenmodes even 
in the dimension(s) of weak confinement will begin to 
change from what they are in effective ID or 2D. 



which has the eigenvalues and corresponding eigenstates 
(unnormalized) 
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VI. CONCLUSION 

We have developed a simple model based on a varia- 
tional approach which can accurately describe the prop- 
erties of a Bose Einstein condensate from the 3D Thomas- 
Fermi regime through various degrees of anisotropy well 



into the regimes of effective ID and 2D. The advantage 
of the model lies in its simplicity, so that analytic ex- 
pressions could be derived for several physical quantities, 
and that it can explain some of the results of experiments 
on such systems without having to resort to complicated 
numerical computations. 

In particular we have found expressions for the chem- 
ical potential which are valid for cylindrical condensates 
for all degrees of anisotropy even where the Thomas 
Fermi expression is completely inadequate. We have ob- 
tained expressions for the total energy and the release 
energy which are valid in 3D as well as in effective ID 
and 2D, and the release energy was shown to agree well 
with experimentally measured values. Our results lead to 
analytic expressions for the condensate wavefunction and 
density, as well as the condensate dimensions for differ- 
ent aspect ratios in cylindrical geometry. We have also 
gauged the variation of the quasi-particle energy spec- 



trum as the condensate deviates from effective ID or 
2D, which should provide a better estimate of spectrum- 
dependent physical quantities in such regimes. Our vari- 
ational ansatz should be useful in studying 2D lattices of 
effective ID condensates such as described in pl|; each 
of those ID condensates also has an axial confinement 
and hence is exactly the type of high aspect ratio cigar- 
shaped condensates we consider here. A study of such 
lattices is a topic of some of our ongoing study. 
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